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1. INTRODUCTION 

Optical interconnects were proposed to relieve the predicted bottleneck in many high-speed digital 
computing systems in the era of big data. This is due to their high capacity, high speed, wide bandwidth, and 
low power consumption [1]—[3]. Different alternatives of optical interconnect technologies can be evaluated 
using the capacity which is directly related to the bandwidth density. The interconnect capacity is defined as 
the time space density product evaluated at a fixed operating power [4]-[7]. A significantly large interconnect 
capacity can be achieved by using a number of spatial channels of free space optical interconnects (FSOIs) 
[8], [9]. FSOIs using vertical cavity service emitting lasers (VCSEL) have witnessed exponential growth in 
vital applications, including data communications, optical sensing, laser printing, and illumination [10]—[14]. 

Like any wireless communication channel, FSOI channels suffer from noise [15]—[17]. In this paper, 
we consider two main sources of noise which are the crosstalk among the channels and the detector amplifier 
thermal noise. Crosstalk arises from the laser beam divergence, and it increases with the density of the 
channels. On the other hand, thermal noise, which is characterized by a noise equivalent power (NEP), arises 
due to the large bandwidth of the FSOI channels [4], [18]—[21]. These two types of noise directly affect the 
speed and bandwidth density of optical channels; hence, the use of error correction coding schemes is 
proposed to combat the noise effects. Several researchers have investigated the use of error correction codes 
(ECC) to optimize the time and space resources of FSOI systems. Neifeld and Kostuk [18] analyzed the use 
of Reed Solomon codes to optimize the channel spacing and data rates of optical interconnects. The results 
were demonstrated for both rectangular aperture and Gaussian propagation models. However, the authors 
overestimated the crosstalk noise power in their work which led to suboptimum findings. Al-Ababneh and 
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Tawalbeh [22] studied the effect of fill factor of the photodetector array and the code rate of the ECCs on the 
bandwidth density of FSOI channels. It presented analytical and simulation results that demonstrate the 
importance of optimizing the fill factor which is defined as the ratio of the used detector diameter over the 
interconnect spacing. This ratio significantly affects the noise and signal powers, so the bandwidth density. 
The implementation of extended Golay codes to model forward error correction in the operation of optical 
data links is investigated in [23]. The authors used Golay codes to compensate for a single and multiple 
inoperative sublinks as erasures. Jun et al. [21] made a comparison among a variety of error control schemes, 
including error detection (ED) and error correction (EC) schemes, that can be used in optical interconnects. 
The comparison was made in terms of area, latency, and power consumption. In this paper, we study the 
performance of FSOIs system model in the presence of, besides the fundamental mode, higher order 
transverse modes. In particular, the effect of interconnect spacing, signal to noise ratio (SNR), and fill factor 
(FF) on the bandwidth density of the FSOIs system is presented. In addition, the impact of using Reed- 
Solomon ECCs (RS-ECCs) of several codeword lengths on the bandwidth density and misalignment between 
the transmitter and receiver arrays is considered. We use the RS-ECCs at the optimum code rate of every 
used codeword length. The rest of the paper is organized as follows: in section 2, we provide a description of 
the FSOIs system and derivations of the bandwidth density of the fundamental and the higher order 
transverse modes. The performance of the studied transverse modes with respect to channel spacing, SNR, 
misalignment, and coding is provided in section 3. Finally, the conclusions and future works are drawn in 
section 4. 


2. RESEARCH METHODS 

For the purpose of the study presented in this paper, we consider an array of nine elements of 
vertical cavity surface emitting lasers (VCSELs) as a transmitter and an array of nine photodetectors as a 
receiver. The interconnect distance between the arrays is set to be D in Figure l(a). The VCSELs are 
arranged as shown in Figure 1(b), where A is the VCSELs’ diameter and B, shown in Figure 1(c), is the 
photodetectors’ diameter. The A parameter represents the distance between any two consecutive 
photodetectors (PDs) and is called channel spacing, where the ratio between B and A is defined as the fill 
factor (FF) parameter. 
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Figure 1. Two-dimensional schematic representation of (a) the FSOIs system model, (b) VCSELs array, and 
(c) photodetectors array 
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In this paper, we have assumed that the parameters of all VCSELs are the same. Moreover, we 
consider only the two major noise sources that are the optical crosstalk noise and detector amplifier thermal 
noise. The optical crosstalk noise results from the divergence nature of the laser beam as illustrated in 
Figures 1(a) and 1(c), while the detector amplifier thermal noise which is also known as the noise equivalent 
power (NEP) characterizes the employed photodetectors. As stated in [22], the bandwidth density (BD) of 
FSOI channels depends on the ratio of the data rate (R) of a single optical interconnect and the channel 
spacing (A) as (1): 


R 
BD= =, (1) 
where R is chosen such that the transmitter (modulator) can achieve this rate and the receiver (detector) can 
demodulate the laser beam efficiently. The relationship between the BD and signal to noise ratio (SNR) as 
well as other parameters of FSOIs system can be explored from the definition of R, which is defined as (20) 


[22], 


(Po XS—Po XN XSNR)? 


R= 
SNR2XNEP2 


À (2) 
where P, is the transmitted power per optical channel in watt, S is the useful received optical signal by 
intended receiver, N is the optical crosstalk noise, and NEP is the noise equivalent power in watt per square 


root of Hertz. The definition of R in the above equation is derived based on the definition of SNR given 
by (3) [22]. 


SNR = 10 log (z) (3) 


PoXN+NEPXR?/2 


It is clear from (1) and (2) that the BD can be found by finding the values of the optical signal power 
(S) and the optical crosstalk noise (N) for predetermined values of P,, NEP, SNR, and A. The value of S 
received by the photodetector PD9 in Figure 1(c) from the VCSEL9 in Figure 1(b) is determined by 
integrating the two-dimensional normalized optical field over the surface of PD9 which is located at plane 
Z = D from the VCSELs array. If we assume that the VCSEL9 is operating in the transverse electromagnetic 
mode mn (TEM), where m and n are integers that represent the radial and angular mode orders, 
respectively, then S can be determined as (4). 


A 
Smn = Sop, Imn (1,9, D)rdrdð = S SZ Imn (r, 0, D) rdrdé. (4) 


With the help of [19], the normalized optical field (J) of a transverse mode TEMmn at a distance D 
from the transmitter is defined as (5), 


hep) Ts (n H) x exp (Fe 


mxWZx(m+n)!x(Sont+D Wp 


) x cos*(nd), (5) 


where Wp is the laser beam spot size at the photodetector that is away distance D from the VCSELs array, 
and is given by: 


2mxAxD\2 1 
A2 , 


Wp =W; ( 


where W; and A are the beam waist at the output of the VCSEL which is equal to A/(V27) and the operating 
wavelength, respectively. The Son function is defined as: 


1, n=0 
Son = P otherwise’ 
and L? (æ) is a Laguerre polynomial of order (m, n) and is given by [24], 
-n m 


a d 
I" (a) = ar x exp(a) x aa (a™*" x exp(—a)). 
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In this paper, we consider three modes of operation for the laser beam which are the TEMoo, TEM 10, 
and TEM». Moreover, we assume that there is no coupling among the laser modes. The normalized optical 
field (5) of the aforementioned transverse modes, respectively, is given as (6) to (8), 


2 

-2r2 2 

ae) © 
2 2 À 

_ 2 —2r o [2r 

= sp (og) eG) o 


aee) o 


where L3 (æ) = 1, L? (a) = 1 — a, and L9 (a) = 1 — 2a + (1/2)a?. 

Upon substituting the Im,n of (6), (7), and (8) into (4), we find the optical signal power (S) received 
by the intended photodetector, PD9, due to the operation of each mode separately. The optical crosstalk noise 
(N) which results from the divergence of the laser beam, as shown in Figure 1(a) and Figure 1(c), can be 
calculated based on the power falling from the VCSEL9 on the photodetectors PD1, PD2, PD3, PD4, PDS, 
PD6, PD7, and PD8. For the fundamental order Gaussian mode, TEMoo, the Nog can be defined as (9), 


2 27? 
Noo = 7A J exp A rdrd@ = 


-H f exp (2 f) raras + 2h "Few (2 


xw 


loo = —— x ex 
00 xxwe p 


and 


bo = —— x ex 
20 axwZ p 


2) rare, (9) 


where rë = (rcos(@) — A)? + (rsin(@))? and r? = (rcos(@) — A)? + (rsin(@) + A)?. In (9), we exploited 
the fact that the optical crosstalk noise received by the photodetectors PD1, PD3, PD5, and PD7 is the same, 
as well as the one received by the photodetectors PD2, PD4, PD6, and PD8 is the same. This is due to the 
symmetry in the distribution of the photodetectors and the rotational symmetry of the laser beam. 

For the TEM 9 and TEM29 modes, the N of each of these modes can be determine as (10) and (11): 


8 2 
N A > L e) x (: a) drd@ 
= ——__,, exp | —— ——,] rdrdé = 
0 ax We é w? w2 
4x2 p? 2r2\" 
E i lk exp|—-] x| 1- rdrd@ 
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(10) 
and 
2 
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where r and r? are same as the ones that are defined in (9). 

Equations (4), (9), (10), and (11) can be solved numerically or by following analytical models to 
solve integrals including the ones proposed in [16]. Having determined Smn and Nmn of each transverse 
mode, the BD can be obtained for any given SNR level and a predetermined values of P, and NEP. The BD 
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of the FSOIs system and the misalignment tolerance between the transmitter (VCSELs) and the receiver 
(PDs) arrays can be significantly improved with coding. For this purpose, we introduce the use of error 
correction codes (ECCs) in the considered FSOIs system. In particular, we use Reed-Solomon (RS) ECCs of 
various codeword lengths to demonstrate the effect of coding on the FSOIs system [25]. The effect mainly 
reduces the SNR level required to achieve a predetermined bit error rate (BER) of the system. The BER is 
defined as in (12) [21], [22], [26], [27]: 


BER = (1 = ert(2 =) (12) 
2 2 2 


where erf in (12) is the error function defined as erf(x) = = i exp (—y?)dy. 


The reduction in the SNR level due to coding is translated into a gain in the bandwidth density, and 
this gain is given by [18], 
_ BDef—BDy 


Ggp = ———., (13) 


BDy 


where BD, and BD, are the uncoded bandwidth density and the effective bandwidth density, respectively. 
The BDey is defined as the code rate (r) multiplied by the coded BD¿, where r is the proportion of useful 
information bits of each generated codeword. 


3. RESULTS AND DISCUSSION 

In this section, we present how design parameters of FSOIs system affect the BD and the 
misalignment tolerance. In particular, we demonstrate the effect of channel spacing, SNR, fill factor (FF), 
and coding on the BD and misalignment tolerance on the system. The design parameters of the simulated 
FSOIs system presented in this paper are set to be as the following: operating wavelength (A)=0.85 um, 
transmitted power per channel (P, )=50 uW, interconnect distance (D)=20 mm, VCSELs’ diameter (A)=0.03 
mm, signal to noise ratio (SNR)=10 dB, and noise equivalent power (NEP)=0.1 nW/,/ Hz. All results shown 
in this section are based on these aforementioned parameters. 

Figure 2 demonstrates the bandwidth density (BD) of the TEMoo, TEM, and TEM29 operating 
modes as a function of channel spacing using FF values of 0.8 (B=0.8A, where B and A are the 
photodetectors’ diameter and channel spacing, respectively), and 1 (B=A). From the results shown in the 
figure, we can see that there is an optimal channel spacing (Aopt) for each FF of each studied TEM mode. For 
TEMoo mode as shown in Figure 2(a), the A,,; that achieves the maximum BD at FF of 0.8 is 0.4 mm and at 
FF of 1 is 0.44 mm. At small values of channel spacing (< 2.6 x 1074m), the BD is very small compared to 
the optimal BD, therefore, it did not appear on the figure. For TEM;0, the values of Aopt at FF values of 0.8 
and 1 are 0.67 mm and 0.65 mm, respectively, as demonstrated in Figure 2(b). Figure 2(c) shows the results 
of TEM operating mode. It is clear from the Figure that Aopt at FF of 1 achieves higher BD than Appt at FF 
of 0.8. In addition, the figure shows that the values of Aopt at FF values of 0.8 and 1 are 0.89 mm and 0.81 
mm, respectively. The dependence of BD on the FF is considered in Figure 3. It is clear from the figure that 
there is an optimum FF value that can achieve a maximum BD for each operating mode. For instance, the 
optimum FF value of TEM; mode is 0.93 which allows the BD of the FSOIs system to be as high as 
5.03 x 10!°bit/sec./m?. 

Figure 4 depicts the BD as a function of SNR levels. The BD of the considered operating modes is 
computed at the optimal value of channel spacing using an FF of 0.8 for each mode, see Figure 2. It is 
obvious from the findings that the BD is highly related to the used SNR level and can be approximated by the 
relation R œ 1/SNR?. This suggests the use of ECCs to improve the performance of FSOIs systems in terms 
of BD. Figure 5(a) demonstrates the coding gain of RS-ECCs of various codeword lengths as a function of 
the code rate. The main idea of exploiting coding gain to improve the BD is revealed by the following 
example. If we consider a codeword of length 127 and a code rate of 0.64, then from Figure 5(a) we find that 
the coding gain is 2.7. This means that the SNR used to generate the results of Figure 2 can be reduced from 
10 to 3.7, 10 + 2.7 = 3.7, while maintaining the same BER. As mentioned in (13), this reduction in the SNR 
level is turned into a BD gain. Figure 5(b) shows the correspondence between the BD gain and the code rate 
for five codeword lengths with SNRyncodeq = 10 dB. It is clear from the results that the coding rate is a 
crucial parameter which directly affects the BD gain due to coding. 
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Figure 2. Bandwidth density versus channel spacing using FF values of 0.8 and 1 for (a) TEM o, (b) TEM 0, 
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Figure 3. Bandwidth density of TEM00, TEM10, and 
TEM20 operating modes versus fill factor 


Figure 4. Bandwidth density as a function of signal 
to noise ratio (SNR) 


The effect of misalignment between the VCSEL array and photodetector array on the BD with and 
without coding using an FF of 0.8 is considered in Figure 6. As expected, the use of RS-ECCs improves the 
BD of the FSOIs system of all operating modes. The improvements vary based on the used codeword length, 
where we note that as the codeword length increased, the improvement in the BD also gets increased. 
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However, the misalignment tolerance at the used FF is the same before and after coding. For example, the 
BD of TEMoo mode Figure 6(a) when the misalignment is zero is identical to the BD when the misalignment 
is 1.14 mm before and after coding. The results of TEM ; and TEM exhibit the same performance as TEMoo 
and are shown in Figures 6(b) and 6(c), respectively. These results are new and have not been reported in the 


literature before. 
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Figure 5. RS-ECC of codeword lengths 15, 31, 63, 127, and 255 effects on (a) coding gain and (b) bandwidth 
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lengths 15, 31, 63, 127, and 255 for (a) TEMoo, (b) TEM 70, and (c) TEM operating modes 
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In order to examine the effect of the FF on the misalignment with and without coding, we generated 
the results shown in Figure 7. One can see from the results that as the FF increases the misalignment 
tolerance to achieve a given interconnect capacity also increases, but this continues until reaching an 
optimum FF. Beyond the optimum FF range, the misalignment tolerance starts decreasing. This behavior is 
common for all examined modes of operation, TEMoo in Figure 7(a), TEM in Figure 7(b), and TEM in 
Figure 7(c). Exploiting RS-ECCs improves the misalignment tolerance of the FSOIs system as demonstrated 
in Figure 7. It is clear from the figure as the codeword length increased, the misalignment tolerance of the 
system over the optimum FF range also increased for some modes. Table 1 summarizes the results of the 
optimum FF ranges and the misalignment tolerances of the studied modes before and after coding of various 
codeword lengths. It is worth mentioning that these results are new and have not been published in the 
literature before. 
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Figure 7. Misalignment tolerance as a function of FF before and after applying RS-ECC of codeword lengths 
15, 31, 63, 127, and 255 for (a) TEMoo, (b) TEM 70, and (c) TEM operating modes 


Table 1. Optimum FF range of FSOIs system for maximum misalignment tolerance in mm 
before and after coding 


Codeword length TEM00 TEM10 TEM20 
(n) Opt. FF Max. misal. Opt. FF Max. misal. Opt. FF Max. misal 
- 0.900-0.925 1.44 0.950-0.956 1.58 0.937-0.963 1.47 
15 0.906-0.938 1.46 0.956-0.969 1.59 0.944-0.969 1.48 
31 0.919-0.938 1.47 0.950-0.969 1.59 0.944-0.975 1.48 
63 0.919-0.950 1.47 0.950-0.975 1.59 0.956-0.963 1.49 
127 0.925-0.938 1.48 0.950-0.975 1.59 0.956-0.969 1.49 
255) 0.925-0.944 1.48 0.950-0.975 1:59 0.956-0.969 1.49 
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4. CONCLUSION 

In this paper, an analytical model for the bandwidth density of a short range FSOIs system with and 
without coding at a given BER is provided. In particular, we demonstrate the effect coding and fill factor on 
the performance of bandwidth density and misalignment tolerance of the FSOIs system. Both the 
fundamental order Gaussian mode of the laser beam and higher order transverse modes have been considered 
in the system design. The results demonstrate that bandwidth density and misalignment tolerance could be 
optimized by selecting appropriate code rate and fill factor levels. This study can be extended to include the 
rectangular propagation model of the laser beam and to consider lenses in the design of the FSOIs system. In 
addition, different types of coding schemes can be considered. 
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